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We exhibit a situation in whicli cosmological perturbations of astrophysical relevance propagating 
through a bounce are affected in a scale-dependent way. Involving only the evolution of a scalar 
field in a closed universe described by general relativity, the model is consistent with causality. 
Such a specific counter-example leads to the conclusion that imposing causality is not sufficient to 
determine the spectrum of perturbations after a bounce provided it is known before. We discuss 
consequences of this result for string motivated scenarios. 
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PACS numbers; 98.80.Cq, 98.70.Vc 

It was recently acknowledged that perhaps the most 
fundamental key questions of string or M-theory could 
be addressed in the framework of time-dependent cosmo- 
logical background Q . Conversely, it was found that this 
new physics may imply new solutions for the dynamics 
of the early Universe. It is the case in particular of the 
pre big-bang paradigm or the more controversial 4] 
ekpyrotic scenario 'J^. In both these last two models, 
the effective four dimensional theory presents a transi- 
tion between a contracting and an expanding phase, i.e. a 
bounce 

In general, both the contracting and expanding phases 
can be described by well controlled low-energy physics for 
which high energy corrections (the a' terms [l| in string 
theory for instance) are negligible. As a consequence, be- 
fore (and after) the bounce takes place, it is possible to 
calculate unambiguously the spectrum of gravitational 
perturbations. On the contrary, the bounce itself de- 
mands full knowledge of such corrections which are, in 
practice, either difficult to implement or simply un- 
known. This is why propagating perturbations through 
the bounce represents a technical challenge but is of ut- 
termost importance, in order to make contact with ob- 
servational cosmology. 

The impossibility to know precisely the dynamics of 
the bouncing phase has led to postulate that its effect 
would reduce at most to a scale independent modifica- 
tion of the overall amplitude 0,0. As a matter of fact, 
both during the radiation to matter [lOl | transition and 
through the preheating which are two other exam- 
ples of short duration cosmological transitions, the spec- 
trum is propagated in a scale-independent way. The pur- 
pose of this letter is to re-examine this assumption in the 
context of bouncing universes. We show, by means of a 
specific counter-example, that such a cosmological tran- 
sition can affect large wavelengths in a scale-dependent 
way and this without violating causality. Thus, this last 
argument cannot be invoked to justify the abovemen- 
tionned postulate. 

The simplest way to study a nonsingular bounce 
is to consider general relativistic Friedmann-Lemaitre- 
Robertson- Walker (FLRW) models with closed spatial 



sections and a scalar field 0- We choose 77 = to be 
the conformal time at which the bounce occurs. Then, 
without loss of generality, the FLRW scale factor in the 
vicinity of the bounce (i.e. for rj ^ 1) can be described 
by a power series expansion in 77 with one free parameter 
for each term of the series. We choose those such that 
the scale factor can be approximated by 
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The first free parameter ao represents the value of the 
scale factor at the bounce. The second free parame- 
ter, namely 770, gives the typical timescale of the bounce 
(the physical duration of the bounce is io — 0.0V0) 
and determines the corresponding de Sitter-like tangent 
model jl2|. the deviation from which being measured by 
the third and fourth free parameters S and f , hence the 
non intuitive coefficients in Eq. they are, in some 
sense similar to the slow-roll parameters of inflationary 
cosmology although they are not restricted to small val- 
ues. It has been shown in Ref. that one should 
restrict attention to 11/5 < C ^ —0.1. If 5 = 0, the 
bounce is symmetric. For |7//7;o| S> 1, the above model is 
assumed to be connected to other cosmological epochs. 
Note also that adding more terms in the series does 
not change the following argument [l^ . 

The free parameter 770 is directly connected to the null 
energy condition (NEC) p + p > 0, where p and p are 
respectively the energy density and the pressure. Indeed, 
at the bounce, Einstein equations imply that 
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where the parameter T is defined by the following ex- 
pression 
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where Lp is the scalar field driving the bounce, behaving 
as Lp — -\- PqT] + ■ ■ ■ (a prime denotes a derivative with 
respect to conformal time) around the bounce, and rUp^ is 
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the Planck mass. Combining both equations, we see that 
the NEC is satisfied at the bounce provided |77o| > 1, 
the limiting case tjq — 1 corresponding to the vacuum 
equation of state p = —p. The previous considerations 
demonstrate that the bounce cannot be made arbitrarily 
short if one wants to preserve the NEC. Under this last 
condition, the short bounce limit is not 7?o ~* but rather 
rjo ^ i, i.e. T ^ 0. Note that since T is related to the 
kinetic energy density of the scalar field to the Planck 
density, one expects T ^ 1. 

Let us now turn to the study of the cosmological per- 
turbations around the previous model (see Ref. [lOj)- The 
gravitational fluctuations are characterized by the curva- 
ture perturbations on zero shear hypersurfaces, i.e. the 
Bardeen potential whose master equation of motion 
reduces to that of a parametric oscillator given by 

u" + [n{n + 2)- Vu (t?)] = 0, (4) 

where the factor n{n+2) arises from the eigenvalue of the 
Laplace-Beltrami operator on the closed spatial sections 
(hence n is an integer). In the above equation, we have 
introduced a new gauge-invariant quantity, u, related to 
$ by (see Ref. for details) $ = V3Knu/{2a'^e), where 
K — 87r/TOp^ and Ti. = a' /a. The quantity 9 is defined 
by 0-^ = a'^{l+p/p){l + n^^). This quantity only 
depends on the scale factor and its derivatives (up to the 
fourth order). 

It is well known that only those modes having n > 1 
are not gauge modes. A crucial point is that the values of 
n of astrophysical interest today are such that n ^ 1 (e.g. 
for finow 1-01, 60 < n < 6 X 10^). They are related 
to the more usual wavenumber k^, commonly used in 
inflationary cosmology, by the relation k"^ = n{n + 2){n — 
1) which shows that, in this last context, fc^ <c 1 (since 

~ 1 at the end of inflation with a very high accuracy) . 
The effective potential Vu (rj) in Eq. can be expressed 
as 

K(r/) = y +3(l-c2), (5) 

where = p' / p' which, in some regimes, can be inter- 
preted as the sound velocity. Contrary to the flat case, 
the effective potential cannot be cast into the form of the 
second derivative of a function over that function. 

The fact that the effective potential only depends on 
the scale factor and its derivatives means that it can be 
calculated for a general bounce given by Eq. (QJ . Far from 
the bounce, but still with ij < rjo, we have Vu{r]) ~ 4. In 
the vicinity of the bounce and in the limit of a short 
bounce satisfying the NEC, one has 

Let us notice that Vuij] — 0) is an extremum of the poten- 
tial only in the case 5 = 0. We see that the height of the 



effective potential diverges in the limit T ^ 0. There- 
fore, in this limit, the bounce will necessarily affect the 
propagation of all the Fourier modes, regardless of the 
values of the parameters ao, S or ^. In other words, the 
spectrum of the fluctuations cannot a priori propagate 
through the bounce without being changed in the NEC- 
preserved short time limit. This conclusion is generic and 
does not depend on the details of the model. 

Is this compatible with causality? Let us first re- 
call that causality relies on the concept of horizon 
which, requiring an integration over time d^{t) = 
a(t) a~^(T)dr, t; being the initial time (which should 
be ti = — oo in the case of nonsingular cosmology), is a 
global quantity. Assuming from the outset an homoge- 
neous background is problematic whenever the implied 
horizon is finite since in this case there exists scales 
larger than d^; somehow, imposing any condition on 
these scales is acausal. In other words, fixes the scale 
limit below which the theory is physically meaningfuU. 
Technically, this means that whenever a homogeneous 
background is assumed, one should proceed as follows: 
solve, mathematically, the perturbation equations for all 
modes A, and then implement causality by restricting 
attention to those satisfying A < d^. To be able to de- 
cide whether this condition is met, one needs to embed 
the local bounce transition into a complete global model 
providing d^. In conclusion, without knowledge of the 
function a{t) for all i, there can be no fundamental prin- 
ciple which would preclude large scales to be spectrally 
affected by a local effective potential (@J. (Note that 
bounces are usually implemented in order to regularize 
the singularity, leading to a geodesically complete uni- 
verse, and hence an infinite horizon.) 

Another possible source of confusion is the often made 
identification of the Hubble scale £^ = a^/a', a local 
quantity, with the horizon d^. Once a given scale A is 
inside the horizon at some time t = to, it remains so for 
any time t > to because the ratio of the horizon to the 
physical scale at time t is 

_ k /■*» dr k dr 

where k is the comoving wavenumber of the scale under 
consideration. The first term is by assumption greater 
than unity and the second one is positive definite, hence 
the statement. By the same token, any scale which is out- 
side the horizon will eventually enter it later (unless a sin- 
gularity develops before). In constrast, as is well-known 
in the inflationary situation, a given scale can either exit 
or enter the Hubble radius, which can be done because 
A > f jj is physical, contrary to A > djj- the case of 
a bouncing universe, the Hubble radius diverging at the 
bounce, all scales are, at some stage, sub-Hubble. For flat 
spatial sections, a super-Hubble scale also means that the 
mode is below the effective potential of the perturbation. 
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This is no longer the case in the curved bouncing model 
for which the mode can be sub-Hubble although potential 
dominated, as was discussed in Ref. 12]. 

We now calculate explicitly the effect of a symmetric 
bounce {S = 0) on the evolution of the cosmological per- 
turbations. In this case, it turns out that the effective 
potential can be expressed as the ratio of two polynomi- 
als of order 24; it is represented in Fig. ^ Obviously, for 



where the constant Cr is given by Cx = [— Stt^^/ {8T)]^^^ 
and where the function Ax(?y) is a representation of the 
Dirac J-function, i.e. limx^o^T('7) = ^iv)- Note that 
even though the potential appears singular, this is noth- 
ing but a computational trick allowing an easy derivation 
of the resulting spectrum (recall that in realistic situa- 
tions, the parameter T must remain small but nonvan- 
ishing The equation of motion 10} of the quantity 

u can now be written as 
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FIG. 1: The effective potential Kj(r;) for the perturbation 
variable w(??) as obtained by using the quartic expansion of the 
scale factor. The values 770 = 1.01, (5 = 0, and ^ = —2/5 have 
been used to derive this plot. The almost undistinguishable 
dot-dashed curve represents a rational approximation to this 
potential which is valid in the vicinity of the bounce and that 
is used in Ref. [T^. Clearly, the potential is well-behaved at 
all relevant times, and so is the corresponding variable u. 

such a complicated potential, the equation of motion Q 
cannot be solved analytically. However, in the case we 
are interested in, namely T ^ 1, the potential has a sim- 
ple behavior, which can be studied simply by retaining 
only the smallest order terms in T. 

In the limit where T goes to zero, the extrcmum of 
the potential Vq diverges while its width, characterized 
by ?7z, shrinks to zero, namely 
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This observation is confirmed by a study of the wings in 
Fig. n whose height and position are found to be 



-§+0(T3/^), (9) 



and exhibit therefore a similar behavior in the NEC- 
preserved short time limit. The previous properties sug- 
gest that we deal with a distributional effective potential, 
and indeed, a more careful analysis reveals that [l^ 
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which is nothing but a Schrodinger-like equation for a dis- 
tributional potential. Right before (superscript <) and 
after (superscript >) the point 77 = 0, but still within 
the bounce epoch, a Fourier mode does not interact with 
the barrier since the potential vanishes and the solutions 
are just linear combinations of plane waves e'v^^*+2)'' 
and e-V"("+2)') with coefficients A> , B> and A<, B< 
before and after the bounce respectively. 

In order to calculate what is the spectrum after the 
bounce, being given some initial conditions before the 
bounce, one must apply junction conditions. In the case 
at hand, the matching conditions are [u] = u(0+) — 
m(0~) = and [u'] — — Cxu(O), the last one coming 
from an integration of the equation of motion in a thin 
shell around 77 = 0. This reduces to 
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B> = A< 
B> = A< 



B<, (12) 
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In the limit T — > 0, the constant Cx diverges and there- 
fore the second term in Eq. (|13|l is dominant. Straight- 
forward algebraic manipulation allows us to determine 
the transfer matrix defined by 



A> 
B> 



A< 
B< 



and we obtain the following expression 
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Some comments are in order at that point. First, one 
sees that, as discussed above, the transfer matrix de- 
pends on the wavenumber. Our calculation permits 
to actually predict accurately what the dependence is: 
cx [n{n + 2)]~^^'^. Moreover, the calculation also predicts 
the dependence of the transfer matrix on the parameter 
^ (except in the limit £ — > for which a different calcu- 
lation must be done [l^l)- A point worth mentioning is 
that the overall amplitude diverges as T — > 0. Since u 
is just a mathematically convenient variable, this is not 
necessarily problematic. Indeed, using the relation be- 
tween u and the Bardeen potential $, and the fact that. 
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at the bounce, Eq. ^ holds, one can show that the spec- 
trum of the Bardeen potential is perfectly finite after the 
bounce, even in the limit T ^ 0. Note also that the fact 
that relevant scales may be larger than the duration of 
the bounce itself does not preclude them to be affected 
by the transition. Finally, it is worth mentioning that 
the above result has been recovered in Ref. using a 
different method, including the numerical factor in the 
overall amplitude. 

In more standard situations, even though the ampli- 
tude of $ might change, the curvature perturbation on 
uniform density hypersurfaces, i.e. the quantity called 
Cbst 01 j does not. Indeed, it satisfies 0| 
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which was been shown ^13^] to hold independently of the 
gravitational field equations. Eq. ((TB|l implies that Cbst 
is conserved under the conditions that there is no en- 
tropy perturbation (Jpnad = 0), the decaying mode of $ 
is neglected, and the scales are super- Hubble (fc <C Ti.). 
Since H = at the bounce, the Hubble radius is larger 
than any relevant scale during a finite interval around 
the bounce, i.e. cosmological scales are not large in the 
super-Hubble sense. Through the bounce, moreover, the 
notion of decaying and growing modes is irrelevant. Two 
conditions out of three being violated, Cbst ^'^ reason 
to be conserved, and hence is not convenient for describ- 
ing a bouncing transition. One should note, addition- 
nally, that since the spectrum of is altered, there is no 
reason why Cbst should not be also spectrally distorted 
through a bounce. 

In this letter, we have demonstrated that there is no 
reason to believe that the spectrum of large scale cos- 
mological fluctuations is not affected by a short duration 
bounce, although this is not necessarily the case (one 
could choose for instance 770 ^ 1, however hard this is 
to reconcile with the field theoretical treatment |0| or 
have a "slow-roll" kind of bounce with 770 = 1 Q)- We 
have shown that this occurs when one approaches the 
NEC violation and we have also demonstrated that this 
effect does not violate causality. This result may find im- 
portant applications: although the calculation discussed 
here is based on general relativity, there is no reason why 
causality should act differently in the framework of, say, 
string theory. In string motivated cosmological scenar- 
ios (as for instance in the pre-big bang paradigm 3] or 
in the ekpyrotic case j^), the calculation of the power 
spectrum of cosmological fluctuations is done in the con- 
tracting phase and the predictions relevant for observa- 
tional purposes then stems from the assumption that for 
sufficiently large scales, perturbations are essentially not 
affected by the bounce. Our result indicates that this 
assumption is far from trivial and may challenge the con- 



clusions reached so far in the literature. 
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